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I. INTRODUCTION 

Electronic transport in mesoscopics has been widely 
studied for the last two decades. Particularly due to the 
interplay between the electronic charge and spin many 
novel features have been found, which lead to new inven- 
tions as giant magnetoresistance, tunneling magnetore- 
sistance or the Datta-Das transistor (see e.g. Rfs.ll||2j|3|y 
for reviews). This has opened the field of spintronics with 
possible future applications in the field of quantum com- 
puting. The fact that spin relaxation and dephasing rates 
are rather long compared to charge excitations life-times 
has recently stimulated the investigation of spin trans- 
port in purely magnetic systems without an accompany- 
ing charge transport. Experimentally, anomalous heat 
transport has been reported in (Sr,Ca, La)\/JJui^O\\- 
materials^, which is explained by the specific dynamics of 
spin excitations in quasi-onedimensional systems^. These 
experiments are still in the incoherent regime where the 
mean free path of the spinons is much smaller than the 
system size. In contrast, Meier and Loss£ investigated 
coherent spin transport through mesoscopic ferromag- 
netic and antiferromagnetic Heisenberg spin chains. For 
the case of the isotropic antiferromagnetic spin-1/2 chain 
they found that the spin conductance is quantized in 
units of order (g /is) 2 /h. Matvee - ^ studied spinon trans- 
port through a one-dimensional Wigner crystal, which 
effectively can be described by an antiferromagnetically 
ordered spin chain. Here, anomalous conductance quan- 
tization phenomena occur which might explain experi- 
ments through quantum point contacts^. 
Motivated by these interesting recent results for spin 
transport through one-dimensional spin chains, we ana- 
lyze in this work the magnetic analog of electronic trans- 
port through zero-dimensional quantum dots in the weak 
coupling regime (i.e. weak exchange coupling between the 
magnetic quantum dot and the magnetization reservoirs). 
In particular, we study the spin current through a quan- 
tum dot consisting of two interacting spin-1/2, coupled 
to two ferromagnetic insulators, which are for simplic- 
ity described by a ferromagnetic Heisenberg model (the 
basic transport features revealed here are expected to 



hold for any reservoir involving magnon- like excitations). 
For the same system Wang et ali^ derived a Landauer- 
Biittiker-type formula for spin current transport. How- 
ever, they studied the magnon transport using a magnon 
representation for the dot states and treated the intra-dot 
magnon-magnon interaction in a mean-field approxima- 
tion. In contrast, we study the case of arbitrarily strong 
intra-dot interaction and, using a real-time formalism de- 
veloped in Ref.llll perform only a perturbative expansion 
in the exchange coupling between the quantum dot and 
the magnetization reservoirs. This is the regime where 
in electronic transport interesting effects like Coulomb 
blockade phenomena occur. 

Our main focus lies on the investigation of possibilities to 
switch the spin current by tuning system parameters as 
e.g. an intra-dot magnetic field or the exchange coupling 
between the two spin-1/2. In Rfs. 1 1211 131 it was shown 
that the latter may be modulated, e.g. by application of 
gates. In contrast to electronic transport, we find that 
the specific microscopic coupling configuration and the 
magnon dispersion relation in the reservoirs has a strong 
influence on the spin current. Interestingly, the spin cur- 
rent can not only be dominated by 1-magnon transport 
processes, but also by 2-magnon cycles, leading to specific 
resonances or antiresonances in the spin current as func- 
tion of the intra-dot exchange coupling. Furthermore, we 
discuss two different mechanisms of driving the spin cur- 
rent, firstly a difference of left and right magnon chemical 
potentials, and secondly a magnetic field gradient. We 
find that they are not equivalent but yield different re- 
sults, unlike the counterpart of electrochemical potentials 
in electronic transport. 

The paper is organized as follows. In section [H] we in- 
troduce the model and the extension of the real-time for- 
malism to bosonic transport. Assuming the reservoirs 
to be magnetized in the same direction, subsequently, in 
section IIII Al the spin current is discussed for the case 
of differing chemical potentials, in section 1111 HI for the 
case of a magnetic field gradient applied, in each case for 
two different dot-reservoir coupling configurations. In 
section llll Cl we deal with the case of antiparallel magne- 
tized reservoirs and finally close with some continuative 
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remarks concerning extensions to this work. 



II. MODEL 

Our model consists of a quantum dot, coupled via lo- 
cal exchange to two magnetization reservoirs. The dot is 
made up of two coupled spin- 1/2, whereas the reservoirs 
are assumed to be ferromagnetic insulators, which can 
be described by a ferromagnetic Heisenberg model, see 
Fig.H 

The total Hamiltonian TL for the system can be written as 
sum of the dot Hamiltonian, the reservoir contributions 
and the exchange part, describing the coupling between 
the dot and the reservoirs (in the following called tunnel- 
ing part corresponding to tunneling of magnons instead 
of electrons in the electronic analog): 



T~t — T~tdot + T~Cl + 7~Lr + Ti-T ■ 
The dot Hamiltonian TLdot reads 



H-dot — —Jdot SlS 2 + Bzee 



(1) 



(2) 



where Jdot is the exchange coupling between the two 
spin-1/2, s z = sf + s| and Bz ee represents a homoge- 
neous magnetic field applied on the dotii^ The reservoir 
contributions are given by a Heisenberg Hamiltonian: 



TL^ — J 1 2 ^ S^-jSf 

(hi) 



B r 2J s z 



+ hS? 1 , (3) 



with r £ {L, R}, J > and (i, j) denoting the sum over 
neighboring sites. 1 -. B r represents a homogeneous mag- 
netic field applied on the reservoir r, whereas b\ models 
a remnant of the dot magnetic field Bz ee - For simplicity 
we assume that this remnant field acts exclusively on the 
reservoir sites neighboring the dot Ml The tunneling part 
is given by 



Si S 7 



(4) 



r=L,Ii. i=l,2 



where J % r are the real exchange couplings between the two 
dot spins and the adjacent spins of the reservoirs, in the 
following referred to as tunnel couplings. 
For the case of parallel magnetized reservoirs we as- 
sume the reservoir ground states to be \Sf = —S, Vi) 
(the modifications for the case of antiparallcl magnetized 
reservoirs will be stated in section UlI Cll . Making use of 
the Holstein-Primakoff transformation, TL r may be diag- 
onalized by introducing the magnon creation (annihila- 
tion) operators a\ k (a r k)- 

TL r = u>rk a\. k a r k, u rk — 2JS[l-cos(ka)] + B r , (5) 



with the lattice constant a. Here and throughout the 
paper we set h = &B = !• Furthermore, in the following 



all energies as also the spin current are normalized to JS. 
Rewriting TLt in terms of the reservoir magnon creation 
and annihilation operators yields 



U T = E ( J rk 4 a rk + h. c.) 



(G) 



r.i.k 



with the fc-dependent coupling J % rk oc — J*(l|fc), where 
(1 1 k) represents the amplitude of the magnon wavefunc- 
tion at the dot neighboring reservoir sites. In JBJ we 
have split off the sf 5^-term arising in 10} and have in- 
cluded its ground state average (i.e. replacing S*i — » — S) 
into the dot Hamiltonian. This leads to an effective field 
Hi = J2 r Jr $ f° r the dot states: 

Hdot = -Jdot sis 2 + B Zcc s z + (Hi s\ + H 2 s z 2 ) . (7) 

Thus our final Hamiltonian TL reads 



TL = TLq + TLt , TLq — TLl + TLr + TLdo 



(8) 



The triplet states |T + ) = | T T ) an( l |T_) = |||) are eigen- 
states of TLdot- Within the sub-space spanned by |S) = 
^(IU)-|4T)) and |T > - ^(ITD + IIT)), i.e. the 
states with s z quantum number m = 0, TLdot corresponds 
to the matrix [ J *g AH = (H 1 - H 2 )/2 measures 

the inhomogeneity of the effective field Hi, which origi- 
nates from the tunnel couplings. 

We consider two different mechanisms driving the spin 
current through the dot, which we denote as A/z- and 
AB -configuration in the following. 

For the A/i- configuration the magnetization of the reser- 
voirs, given by the magnon number, is supposed to be 
conserved. This means that we assume the spin relax- 
ation time in the reservoirs to be longer than the typical 
spin current measurement time. Therefore a chemical 
potential \x r is assigned to the magnons in each of the 
reservoirs. Moreover we set B r = here, so that only 
the difference between the magnon chemical potentials 
jil and iir drives the spin current. 

For the AB -configuration we set ji r = in both reser- 
voirs. Here a magnetic field gradient B^/^t) = Bq ± 
AB , switched on simultaneously with the tunneling, 
gives rise to a non- vanishing spin current. The offset 
field Bq > \AB\ introduced here is required for the sta- 
bilization of the ferromagnetic ground state. 

We use the real-time formalism described in Ref. ^] 
to compute the stationary spin current. By integrating 
out the reservoir degrees of freedom this formalism yields 
an effective description of the system in terms of the dot 
degrees of freedom. It is based on the formally exact 
kinetic equations (tunneling and magnetic field gradient 
are assumed to be switched on at time t = 0) 



I r (t) = Tr 



dot 



dt! E Ir (t-t') P (t') 



(9) 



p(t)=-iL oP (t)+ / dt' ^(t-t')p(t'), (10) 
Jo 
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where I r {t) denotes the time-dependent expectation 
value of the spin current through the dot. The current 
operator I r is given by 



I r = -i[H,N r ]_ =i^( J r 



a lk s i 



Jrk4 a rk) (11) 



with [•,•]_ denoting the commutator and N r — 

J2k a tk ark ' so ^ na ^ l ne s P m curren t Ir is positive, when 
a 'spin-up' leaves the reservoir r. p(t) is the time- 
dependent reduced density matrix, which is the trace 
over the reservoir degrees of freedom, Tr res , of the full 
density matrix. In contrast, Tr^ ot denotes the trace over 
the local (dot) degrees of freedom. The free propagation 
of the system is determined by the Liouvillian Lo . It is 
defined by Lo = [Ho, ■ ]_ . The coupling to the reservoirs 
is described by the integral kernels E and E/ r . They are 
given by 

E(i-t') 

= (-i) 2 Tr res [L T e- iLot Te- i ^' ciTLr(T) e iLot 'L T p^ J 



irrcd. 

(12) 



(13) 



irrcd. 



where we assumed that the reservoirs are in equilibrium 
described by the reservoir density matrix p™L- Lt = 
[Ht, ■ ]- is the interaction part of the Liouvillian, and the 
interaction picture is defined by Ly(i) = e lL °'LTe~ lLot . 
The superoperator Aj r is given by the anticommutator 
with the current operator I r : 



A, 



(14) 



T denotes the time ordering operator, and the in- 
dex u irred." indicates that only irreducible diagrams 
are taken into account (meaning that each vertical cut 
through a diagram will cut at least one bosonic reservoir 
line, see Ref. [llj for details). Introducing the Laplace 
transforms f(z) = dte izt f(t) of the time-dependent 
functions f(t) and using the identity lim^oo f(t) 

i lim 2 ^;o+ zf(z) leads to the following expression for the 
stationary spin current If: 



I s 



Tr 



dot 



£/,(io> s 



(15) 



where the stationary reduced density matrix p st is deter- 
mined by 



-iL n 



E(iQ+) 



p st = 0. 



(16) 



We calculate the kernels E and Ej r to lowest order in 
tunneling, i.e. second order in the coupling constants J % r . 
The diagrams contributing are listed in the appendix. A 
typical diagram is shown in Fig. ^a) , which corresponds 
to the following expression: 



(«) 
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FIG. 1: (a) A typical diagram contributing to the kernels E 
and Ej r . (6) The hatched block represents arbitrary contrac- 
tions. 



^E( s il J v«^l s 2>(4l4 /HS 7| s i> 

r(w) n L/R .{u) 



2n 



h3 



(17) 



duj- 



iO' 



where A s ' iS = E„i — E s ^ and E s are the eigenvalues of 
Hdot- n r (u) = 1/ (eP^~ ^ — l) is the Bose function and 
the spectral function T{uj) is given by 



r(w) =%SY)^ \k)\ 2 8{u-u k ), 

k 



(18) 



with u k = 2J5[1 - cos(fca)] + B . In (JHJ TAB corre- 
sponds to the left/right reservoir. Due to the conserva- 
tion of the total spin the following equation holds for the 
s z quantum numbers of all non-vanishing diagrams, as 
indicated in Fig. ^b): 



Ills 



(19) 



Assuming the dot to be in thermal equilibrium initially, 
the non-diagonal matrix elements p(0) s ^ s > of the initial 
reduced density matrix vanish. From Ijl9|l it then follows 
that the only non-trivial non-diagonal elements of the re- 
duced density matrix p(t) (particularly p st ) generated are 
those between the two states with m = 0. 
In some cases the influence of these non-diagonal ele- 
ments onto the stationary spin current may be neglected. 
Then l|15f) and l|16f) reduce to classical master equations: 



r,s' 9=± 



(20) 
(21) 



where q = ± corresponds to transitions with a spin-flip 
of Am = ±1 in the dot. The rates E^, are given by 



v- L /«+ 



~\L/ Ft- 



('\EiJi,*4W)\ 

x r(A ss / =p AB) n L/Ii (A ss > =p AB), (22a) 

H<s'IE^ /R s+| S }| 2 

x T(A S , S t AB) [1 + n L/R {A s , s =f AB)] , 

(22b) 



corresponding to absorption and emission of magnons by 
the dot, respectively. This is the general form for arbi- 
trary AB and fi r (the latter enters the Bose function). 
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FIG. 2: The two coupling configurations studied: (a) Serial 
coupling, (b) Parallel coupling, note the minus sign. 



For investigation of the two driving mechanisms we set 
either AB = (A/i-configuration) or fi^ = fin = 
(A_B-configuration). The product of the spectral- and 
Bose function, entering l|22a|) . is depicted in Fig. [21 (left) . 
While the spectral function introduces the band edges, 
the Bose function leads to the exponential decay of the 
rates. 



III. RESULTS 



To ■ 



Jdot 
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FIG. 3: (a) The dot level structure for |S)/|To)-eigenstates 
of Tidot and its dependence on the dot magnetic field Bzcc 
and the intra-dot exchange coupling Jdot- The average H of 
the effective field Hi adds to the homogeneous magnetic field 
Bzcc- (b) The possible transitions in case of serial coupling and 
|S)/|To)-eigenstates. Each transition indicated by an arrow 
can occur over both of the reservoirs. With increasing Ast_ 
beyond the band edge lol/r > the ST_ -transition over the 
reservoir L/R becomes possible. 



FIG. 4: The largest contributions to the spin current II, 
exemplarily for serial coupling and two different level config- 
urations. Darker arrows correspond to larger currents. 



First we state some basic principles determining the 
single magnon transport through the dot, for the case of 
parallel magnetized reservoirs (the modifications for the 
case of antiparallel magnetized reservoirs will be stated 
in section ITf I C|) . 

Tunneling of a magnon out of/into a reservoir involves 
a spin-flip Am = ±1 in the dot. Moreover it requires 
an transition energy above the lower band edges uj^ /r of 
the corresponding reservoirs^ The magnon dispersion 
relation enters the rates via the spectral function T(lo). 
At last the tunnel couplings J l rl which enter the rates as 
prefactors to the spin matrix elements, impose further 
restrictions onto the tunneling transitions. Fig. [21 shows 
the two coupling configurations studied, which we refer 
to as serial and parallel coupling. For instance these may 
be realized by two spin-1/2 quantum dots coupled to the 
reservoirs in series or in parallel, respectively. Particu- 
larly for the case of parallel coupling we have examined 
the influence of a relative phase between the couplings. 
For vanishing phase no qualitatively new features occur 
in comparison to the serial coupling case. But for a phase 
of 7T new features, such as the 2-magnon tunneling, arise. 
Therefore we restrict the discussion of the parallel cou- 
pling to this case, which is depicted in Fig. Elb). 
Depending on the couplings some states may decouple 
from either of the reservoirs. E.g., in the case of parallel 
coupling transitions between |S) and |T_) can occur over 
the left reservoir only, while transition between |Tq) and 
|T_) occur exclusively over the right one; see Fig.[2fb). 



Fig. E2a) shows the dot level structure, exemplarily 
for |S)/|To)-eigenstates, i.e. AH = 0, and its dependence 
on the dot magnetic field Bz ee and the intra-dot ex- 
change coupling Jdot- The average H = (Hi + H 2 )/2 of 
the effective field Hi adds to the homogeneous magnetic 
field Bzee- The possible transitions in the case of serial 
coupling are illustrated in Fig. |3fb) [note that transi- 
tions can only occur between states where the state with 
higher energy has also a higher spin quantum number 
in z— direction, since absorption (emission) of a magnon 
by the dot increases (decreases) energy and spin simul- 
taneously]. In this case, the couplings impose no further 
restrictions, i.e., if the transition energy between the dot 
levels is sufficient, the respective transition can occur over 
both of the reservoirs. 

In view of the spin current the various transitions con- 
tribute differently; see Fig. where only the essential 
current contributions are shown. Due to the exponential 
decay of the magnon occupation with energy, large tran- 
sition energies lead to small current contributions due to 
suppressed magnon absorption processes (note that cur- 
rent can only occur if there is a closed cycle of magnon 
absorption and emission processes). Therefore the T+S- 
contribution is left out. Furthermore the upper levels 
are generally less occupied and thus play also a minor 
role in current transport, as long as other 'channels' are 
available, as is the case in Fig. 0] Since the stationary 
dot spin is conserved it is sufficient to consider the cur- 
rent II over the left reservoir. W.l.o.g. we choose Hl/r, 



Jdot 



FIG. 5: The spin current I L (xl0 4 /JS) for the A/i- 
configuration in serial coupling as a function of the dot mag- 
netic field Bzcc and the intra-dot exchange coupling Jdot, both 
normalized to JS. The lines indicate the onsets of transitions. 
[T = 0.05, (Jl,Jl, Jr,Jr) = (3,0,0, 1) • 10~ 2 , AH = 0.005, H = 
0.01, 6i = 0.1, hl >r = (-0.01, -0.3)] 
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resp. AS such, that the spin current flows from the left to 
the right reservoir. Due to its insignificance the backflow 
from the right to the left is also not depicted in Fig. 0] 



A. Spin current for A^i-configuration 

Serial coupling. For the A/x-configuration the band 
edges are given by u>l = = 0. Therefore transitions 
are possible for arbitrary small transition energies. Fig.0 
shows the spin current for the A/x-configuration in serial 
coupling as a function of the dot magnetic field Bzee and 
the intra-dot exchange coupling Jdot- With each new 
channel opened the spin current increases. For Jdot — 
the current vanishes, since in this case the whole system 
separates into two uncoupled parts (left and right); see 
Fig. El^a). The cross section indicated by the arrow is 
shown in Fig. where the level schemes illustrate the 
main current contributions. Each time the singlet level 
crosses one of the triplet levels, a qualitative change of the 
spin current characteristics occurs. The outstanding cur- 
rent peak is due to the onset of the ST_-transition, open- 
ing a new current channel which involves the strongly oc- 
cupied T_-level. 

The width of the dip at Jdot = is determined by the in- 
homogeneity AH of the effective field Hi\ see upper inset 
in Fig.EI T° understand this, consider the Jdot = eigen- 
states |T I}/ 1 IT) of Hdot, to which a finite Jdot < AH may 
be regarded as a perturbation. The corresponding level 
scheme is depicted in the lower inset. Here, each tran- 
sition may exclusively occur over one of the reservoirs 
only. As a consequence the stationary current vanishes 
for Jdot = 0, since there exists no path leading back to 
an arbitrary initial state while at the same time carry- 
ing spin from the left to the right reservoir. E.g., on the 
path — > |T + ) — > HI) a magnon tunnels from left 



FIG. 6: The spin current II as a function of Jdot (cross section 
in |SJ . The level schemes at the bottom illustrate the main 
current contributions, compare Fig. 2] Inset: dip for different 
inhomogeneities AH and level scheme in the neighborhood of 
Jdot = 0. [T = 0.05, B Zcc = 0.1, {JI,JI,Jr,Jr) = (3,0,0,1) • 
10~ 2 , AH = 0.005, H = 0.01, bi = 0.1, (j.l,r = (-0.01, -0.3)] 



to right. But in order to get back to the initial state, 
which is necessary to get a stationary current, a magnon 
must tunnel back again. So effectively the stationary spin 
current vanishes. Switching on Jdot evokes a 'coupling' 
between the |||)/|||)-states. Since this coupling opens a 
new channel between |||) and |||) without magnon flow, 
the paths indicated in the inset level scheme give rise to 
a nonvanishing spin current. The J^ot-mediated coupling 
becomes important for Jdot of the order of AH, which 
measures the splitting between the two levels. Thus AH 
determines the width of the current dip. 
For | AH" | <C |min{J*}|, particularly AH = 0, a narrow 
dip still remains. In this case the eigenstates are S/To- 
like also in the neighborhood of Jdot — 0. Then, the 
suppression of the spin current in the dip is due to a 
destructive interference between equivalent current car- 
rying paths, e.g. the ST_- and ToT_-paths. Thus, here 
the nondiagonal matrixelement p|' T has to be taken into 
account, i.e. only the full kinetic equation yields correct 
results in the neighborhood of Jdot = 0. 
Tuning further up to Jdot > Bzcc, the spin current fi- 
nally saturates in the contribution of the ToT_-transition, 
while the S-level occupation gets negligible. 

Parallel coupling. Fig. [7| shows the spin current for 
parallel coupling as a function of Bz ee and Jdot- For 
Jdot = instead of a dip, now a peak occurs in the spin 
current. In contrast to the current peak in Fig. this 
peak cannot be ascribed to the onset of a transition. To 
understand this we discuss the current along the cross- 
section indicated; see Fig. Again, the level schemes 
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FIG. 7: The spin current I L (x itf/JS) for the A/i- 
configuration in parallel coupling as a function of Bzee and 
Jdot- In contrast to serial coupling here a peak occurs for 
Jdot = (note the logarithmic scale). 

[T = 0.05, (Jl, Jl, Jjt, J%) = (1, -1, 2, 2) • 10~ 2 , AH = 0.005, H = 
0.01, bi = 0.1, ii hlR = (-0.01, -0.3)] 
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FIG. 8: The spin current II as a function of Jdot (cross 
section in 0. Note that here AH plays the role of the 
perturbation. The inset shows the magnified sector (a). 
[T = 0.05, Bzco = 0.15, {Jl,Jl, Jr, Jr) = (1, -1, 2, 2) • 10~ 2 , AH = 
0.005, H = 0.01, 6i = 0.1, Mi. J? = (-0.01, -0.3)] 



illustrate the largest contributions to the spin current. 
The difference between serial and parallel coupling is, 
that in the case of the latter the transitions are more re- 
stricted [as mentioned in the context of Fig. 0[b)] ■ Tran- 
sitions involving |S) (|Tq)) may occur exclusively over the 
left (right) reservoir here. Therefore a stationary spin 
current requires a sequential tunneling of two magnons, 
e.g. on the path |T_) — > |S) — » |T + ) — > |T ) — > |T_), as 
illustrated in the level scheme for sector (d). 



The current peak (c) appears for Jdot since then 

the inhomogeneity AH becomes essential by 'coupling' 
the |S)/|To)-states and thus introducing a further channel 
which enables 1-magnon tunneling; sec the central level 
scheme. In this way the otherwise inevitable 2-magnon 
tunneling is cut short, resulting in an increased current 
[compared to the 1-magnon tunneling peak (c) the cur- 
rent in sector (d) is exponentially decreased by a factor 
ofexp(-/?A To T_)]. 

Unlike in sector (d) the AH-mediated coupling still takes 
influence in sector (b) , as illustrated in the according level 
scheme. While in the former case the ST_ -transition is 
the 'bottleneck', which cannot be short-circuited, in the 
latter case the restricting T+S-transition is cut short in- 
deed. 

Besides these qualitative arguments it shall be mentioned 
that the quantitative computation of the spin current is 
performed by setting up the kinetic equation in terms of 
the eigenbasis of TLdot ■ The influence of nondiagonal ma- 
trixelements between the two states with m = 0, which 
gets most important for a small level splitting, has to be 
considered too. For Jdot — the eigenstates are |TI)/|lt) 
while their minimal splitting is given by AH. Since the 
inhomogeneity AH cannot vanish for the J* in parallel 
coupling, the states are never degenerate and thus a de- 
structive interference is largely suppressed. However a 
small current correction still remains, resulting in a re- 
duced peak height. Thus, in order to account for this 
correction, the full kinetic equation has to be considered 
instead of the classical master equation. 



B. Spin current for AB-configuration 

In the case of the AB-configuration a magnetic field 
gradient drives the spin current. Here the minimal exci- 
tation energies in the left and right reservoir differ, since 
the band edges are shifted by the magnetic field gradi- 
ent: Ul/r = Bo ± AB. Therefore, in contrast to the 
AyU-configuration, the transitions over the left and right 
reservoir set in for different dot transition energies which 
leads to additional features in the spin current. The en- 
ergy dependence of the rates entering the kinetic equa- 
tion, particularly the shift of the band edges, is shown 
in Fig. It also illustrates the successive onset of the 
transitions over the right and left reservoir, exemplarily 
for the ST„ -transition. 

Serial coupling. Fig.llOlshows the spin current for the 
AB-configuration in serial coupling. Since we choose 
AB = Bq here, it follows that lor — 0. However, a 
different choice does not lead to relevant modifications. 
Again the onsets of transitions are indicated by lines. 
For Jdot S Aff deviations from these lines occur due to 
the influence of the inhomogeneity Aff , which shifts the 
eigenenergies of TLdot non-linearly. 

In sector (I) the current is basically the same as for the 
case of the A/x-configuration, compare Fig. [S] Beyond 
this similarity new features arise, specific to the AB- 
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FIG. 9: (Left) The product of spectral- and Bose func- 
tion, which determines the energy dependence of the rates 
[compare JHJ], f° r the left and right reservoir as a func- 
tion of transition energy Asj_ . The band edges are given 
by u-> L / R = Bo ± AB. (Right) The onset of transitions illus- 
trated by means of the ST_ -transition; for clarity the upper 
levels are left out. When Ast_ exceeds ujr/l the transition 
over the right/left reservoir becomes possible. 




FIG. 10: The spin current I L (xlQ 5 / JS) for the AB- 
configuration in serial coupling as a function of Bzcc and Jdot- 
[T = 0.05, (Jl, Jl, J R , J%) = (2,0,0, 1) • 10~ 2 , AB = 2.5 10" 3 , H = 
7.5 • 10~ 3 , bi = 0.1, AB = B = 0.05] 



configuration. In the lower half of Fig. ^| an area of zero 
current occurs, denoted as sector (II). The left reservoir 
is completely decoupled here since all transitions involv- 
ing it are disabled (all the transition energies lie below 
the left band edge). Furthermore there are additional 
steps in the sectors (III) and (IV). To understand these, 
consider Fig.llllwhich shows the cross section indicated. 
In sector (III) the onset of the ST_-transition over the 
right reservoir [compare Fig. Efb)] does not open a new 
current channel but leads to a swap of occupation be- 
tween the S- and T_ -level; see inset. Thus the increase of 
the spin current is due to the enhanced T_-occupation. 
This is revealed by the fact that the T_-occupation is ef- 
fectively mapped onto the current (compare inset). An 
analog argumentation holds for the step in sector (IV). 
In contrast to this, the peak (I a) arises since a new (res- 
onant) current channel is opened by the onset of the 
ST_-transition over the left reservoir. This is completely 
analog to the corresponding peak in the case of the A/i- 
configuration (compare Fig. ■ F° r sectors (lb) and (Ic) 




FIG. 11: The spin current II as a function of Jdot (cross 
section in Fig. I1U1 . In sector (III) the ST_ -transition over 
the right reservoir sets in. Although no new current chan- 
nel is opened, the current increases due to a swap of oc- 
cupation between the S- and T_ -level; compare the inset 
showing the occupation of the S- and T_ -level in this sector. 
[T = 0.05, B z „ = 0.17, (Jl, Jl,Jr,Jr) = (2,0,0,1) • 10~ 2 , AH = 



2.5 • 10" 



H = 7.5 • 10" 



bi = 0.1, AB = B 



: 0.05] 



there is no qualitative difference as well since all transi- 
tions are possible over both reservoirs, as is the case for 
the A/i-configuration too. 

Parallel coupling. Analog to the case of the A/i- 
configuration (compare Fig. [S] and 0) for Jdot = the 
dip is replaced by a peak in sector (I). Nevertheless the 
sector (II) of zero current, as also the steps in sectors 
(III) and (IV) remain as for the case of serial coupling. 

Finally, for the case of serial coupling, we discuss the 
spin current as a function of the magnetic field gradient 
AB; see Fig.^J where the dot parameters Bz cc and Jdot 
are fixed. For small AB the current increases linearly. 
Then, due to the shift of the band edges, the transitions 
over the left reservoir set out successively (see inset), each 
resulting in a drop of the current. When the left reservoir 
finally decouples, the current vanishes completely. 
The shape of the current as function of AB depends 
strongly on the dot level structure determining which 
transitions are involved in which order. Fig. 1131 shows 
exemplarily the current as a function of AB and Jdot- 
Once again with each onset of a transition the current 
increases. 



C. Antiparallel magnetized reservoirs 

So far we have considered the reservoirs to be magne- 
tized in the same direction. In the following we briefly 
discuss the case of two antiparallel magnetized reservoirs, 
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FIG. 12: The spin current II as a function of the magnetic 
field gradient AS. Inset: the transition energies and the band 
edges uJl/r = So ± AS. When a band edge (here u>l) exceeds 
a transition energy the corresponding transition sets out. The 
splitting between Ay + T and At t_ is approximately given by 

2AH 2 1 Jdot = 3.1-10 -4 , which cannot be resolved on the scale 
shown. [T = 0.05, J dot = -0.04, B Z€€ = 0.18, j£, J^, j|) = 

(2,0,0,1)-10 -2 , AH = 2.5-10 -3 , H = 7.5-10 -3 , b t = 0.1, B = 0.12] 
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FIG. 13: The spin current I L (xW 5 /JS) for the AS- 
configuration in serial coupling as a function of the mag- 
netic field gradient AS and Jdot- [T = 0.05, B Z cc = 
0.18, (Ji, JI,Jr,Jr) = (2,0,0, 1) • 10~ 2 , Aff = 2.5 • lO" 3 , H = 
7.5 • 10~ 3 , 6i = 0.1, S = 0.12] 



where w.l.o.g. the ground state of the right reservoir is 
flipped upwards. The real-time formalism used is also ca- 
pable to deal with such a system. This is simply achieved 
by exchanging the magnon creation and annihilation op- 
erators for the right reservoir in JSJ, © and i]l lfipi^ 
In this antiparallcl configuration a spin current flows even 
without an additional driving mechanism, such as a mag- 
netic field gradient applied, since the reservoir magneti- 
zations tend to equalize by exciting magnons in both of 
the reservoirs. Fig. 1141 shows exemplarily the spin cur- 




Jdot/(B Ze e + H) 





FIG. 14: The spin current II as a function of Jdot for the 
case of antiparallel magnetized reservoirs in serial coupling. 
When the S-level lies between the T_- and T+ -levels the cur- 
rent is suppressed due to decoupling of the right reservoir. 
In the remaining sectors the current is carried by the trans- 
port cycles depicted in the level schemes. [T = 0.05, B Zcc = 
0.1, (Jl, JI.Jr, J%) = (3,0,0, 1)-10 -2 , AH = 0.01, H = 0.005, 61 = 
0.1, hl,r = (-0.03, -0.03)] 



rent as a function of Jdot for the case of serial coupling. 
In view of the basic transport principles the main change 
is, that transitions over the right reservoir can now only 
occur between states where the state with higher energy 
has the lower spin quantum number in z— direction, since 
the magnons in the right reservoir carry a 'spin-down'. 
The spin current II is negative since 'spin-ups' are enter- 
ing the left reservoir, in order to equalize the magnetiza- 
tions. The current is suppressed, when the singlet level 
lies between the two outer triplet levels, since then the 
right reservoir is completely decoupled. In the remaining 
sectors a non- vanishing spin current is carried by magnon 
tunneling processes which involve all dot levels, as indi- 
cated in the level schemes in Fig. El Since these pro- 
cesses require a magnon to be initially excited in the left 
reservoir (long upward arrow in the level schemes), the 
current decays exponentially with increasing | Jdot \ ■ Fur- 
thermore, due to this the current vanishes with vanishing 
temperature. 

We conclude that also in this antiparallel configuration 
the spin current can be switched by tuning the intra-dot 
parameters Bz ee and Jdot, e-g- by simply evoking a ST_- 
crossing. 
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Conclusions 

In summary, our investigation of the proposed spin 
quantum dot system revealed a rich structure of the spin 
current as a function of the intra-dot parameters B Zee 
and J dot- Due to the involved large current variations, a 
switching of the current is possible in principle. 
For the case of parallel magnetized reservoirs the influ- 
ence of the dot-reservoir couplings on the spin transport 
becomes apparent by opposing the serial and parallel cou- 
pling configurations. The restriction of the transitions in 
the case of parallel coupling leads to drastically differ- 
ent transport characteristics, which originate from the 
sequential 2-magnon tunneling. However, due to the dif- 
fering signs of the dot-reservoir couplings, the parallel 
coupling configuration represents rather an extreme case 
compared to the serial coupling configuration which sug- 
gests itself. It shall be mentioned that small deviations 
from the 'pure' coupling configurations discussed here, do 
not lead to qualitative modifications of the spin current 
results. 

The concept of a magnetochemical potential as sum of 
the chemical potential and the magnetic field, analog 
to the electrochemical potential in the electronic case, 
is deficient. The magnon chemical potential controls 
the magnon occupation, whereas an applied magnetic 
field gradient shifts the magnon band edges while leav- 
ing the overall magnon occupation unchanged. As we 
have shown, both mechanisms, the differing left and right 
magnon chemical potentials as also the magnetic field 
gradient, yield different results for the spin current. This 
is due to the fact that the magnon dispersion relation, 
particularly the band edges, play an essential role in spin 
transport, in contrast to the case of electronic transport, 
which primarily takes place between the electrochemical 
potentials. 

In the case of antiparallel magnetized reservoirs the 
main features of the spin current can be explained by 
considering slightly modified spin transport principles. 
As we have shown, here spin current switching is possible 
too. 

We remark, that in those sectors where the spin 
current carried by single magnon tunneling is expo- 
nentially suppressed, a cotunneling of magnons, analog 
to its electronic pendant, may yield non-negligible 
contributions to the spin current. However, the current 
suppression due to magnon dispersion effects should 
remain untouched. 



We gratefully acknowledge discussions with F. Meier, 
J. Martinek, D. Loss and H. Capellmann. 



APPENDIX 

The diagrams contributing to the kernels E(iO f ) and 
E/ r (i(F) to lowest order in tunneling are shown in Fig. 1151 
Thereby, in matrix representation the kernels are given 

by 

w^vs = E E ^ (a- 1 ) 

r=L,R d=l 

£i r (itt% lSl , S2 4 = (1) + (2) + (5) + (6) . (A.2) 
By defining the auxiliary function 

r H<H 



^{AE) = ±-jdu 



AE-B r -uj + i(L 



(A.3) 



where p E {+,—}, = n r (u>), n~(u) = 1 + n r (uj) 

and Bl/r = ±AB, the expressions corresponding to 
the diagrams with left running reservoir contractions are 
given as follows: 

(i) = E M&t^MVitjM) <, + ( A W . 
(3) = E (*il^«7l«2>tel A + K> 4"( A 4%) > 

(5) = -EE S « ( s i\ J r4\s) (s\J J r sj\s 2 ) o\- (A s , s ) , 

<7> = -EE^ (4l^7l s >< s l J >MK + ( A -i)- 



1,3 



Each of the remaining diagrams (d) with right running 
contraction is obtained by taking the complex conjugate 
and interchanging the primed by unprimed indices (and 
vice versa) in the expression for the left running diagram 
(d-1). E.g.: 



(2) = 



(s^Ks+ls^is^sJls^al + iA^,) 
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The factor p/ie is included in the Zeeman energy Bzee, in 
the following referred to as a magnetic field. This is also 
the case with all following magnetic fields, resp. Zeeman 



energies. 

In this work we assume one-dimensional reservoirs. How- 
ever, since the qualitative shape of the spectral function 
does not depend on the reservoir dimension, apart from 
minor modulations, no qualitative changes occur for reser- 
voirs of other dimensionality. 

The consideration of the remnant field b\ is necessary to 
break isotropy. This is due to a subtlety occurring in the 
calculation of the spectral function T(u>) [cf. I18fl ] for the 
semi-infinite Heisenberg model. As a consequence of this 
&i = represents an unphysical border case. 
Isotropy may also be broken by considering an anisotropic 
Heisenberg model (z-anisotropy). However, as long as a 
ferromagnetic ground state prevails, this does not lead to 
qualitative modifications of the results presented here. 
For the parameters in consideration the upper band edge is 
of minor importance because there the magnon excitation 
is already suppressed exponentially. 

This results in exchanging <-» and reversing the 
sign of lu — > —lo in the resolvent of the expression JA.3II 
[the latter involves a sign reversal of the arguments of the 
spectral- and Bose functions in 122H 1. 



